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Abstract. We consider a class of evolution equations taking place on the edges of a 
finite network and allow for feedback effects between different, possibly non-adjacent 
edges. This generalizes the setting that is common in the literature, where the only 
considered interactions take place at the boundary, i.e., in the nodes of the network. 
We discuss well-posedness of the associated initial value problem as well as contractivity 
and positivity properties of its solutions. Finally, we discuss qualitative properties that 
can be formulated in terms of invariance of linear subspaces of the state space, i.e., of 
symmetries of the associated physical system. Applications to a neurobiological model 
as well as to a system of linear Schrodinger equations on a quantum graph are discussed. 



The mathematical analysis of elliptic operators acting on spaces of functions on networks 
was started by G. Lumer in [iQ]"!!!]- It has been subsequently continued by many authors, 
both in mathematics (in the context of network diffusion problems, see e.g. |27)-|29j-|25j) 
and in physics (leading to the theory of quantum graphs, see e.g. |ll)-|18j-|17j). 

A form of weak nonlocal interactions for evolutionary problems over network-shaped 
structures has already been considered in e.g. Additionally, we are interested in 

discussing systems of strongly coupled evolution equations. Such couplings may correspond 
to the cases of either a phenomenological interaction among parts of the physical system 
(like in a certain neurophysical theory, which we briefly discuss in Section 5.1) or else as a 
form of external control (possibly with the aim of stabilization) . 

More precisely, we want to allow the evolution in a point of the network to depend nonlo- 
cally on those finitely many other points of the network G that have same parametrization 
with respect to the network edges. In other words, we will discuss the strongly coupled 
elliptic operator defined by 



where Uj represents a relevant physical quantity on the j edge of the network. The operator 
A is the gradient of the energy functional E defined by 

/• 1 ^ 



As usual in the context of evolution equations on networks, we also allow for a further, 
weak form of interaction given by a generalized Kirchhoff-type law in the ramification nodes. 
These two forms of interactions between individual linear elements give rise to a well-defined 
system of diffusion or Schrodinger equations. Dwelling on interesting similarities with the 
biological theory of neuronal coupling (cf. Section 5.1), we often call ephaptic and synaptic 
the influences that depend on the behaviour of the process in another edge or in another 
node of the network, respectively. 



1. Introduction 



(1.1) 




i=l 



X e (0,1), j = 1,..., 
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Well-posedness of such diffusion and Schrodinger problems can be proved under quite 
general conditions on the coefhcients (cji). Instead, further qualitative properties strongly 
depend on the coupling coefficients that are actually considered. In particular, we can show 
that in spite of the parabolic nature of the diffusion problem, no maximum principle holds as 
soon as the ephaptic coupling is nontrivial - i. e., as soon as the matrix {cji) is non-diagonal. 

In the second part of this paper, we discuss the issue of symmetry properties for both 
diffusion and Schrodinger equations on networks. 

One says that a given physical system exhibits a symmetry if some of its properties 
remain invariant under the action of a certain class of transformations. More precisely, in 
the Lagrangian formulation of field theory, one says that there exists a (global) symmetry 
of a given dynamical system if the Lagrangian C{4') of the field 4> is invariant under all 
(time- and space-independent) transformations O that belong to a group O, the so-called 
gauge group of the system, i.e., if C{(f>) = C{04>). The prototypical example is given by 
the invariance under rotations of the Laplacian: this implies a symmetry for both the heat 
and the Schrodinger equations in R", whose gauge group is the orthogonal group O = On- 
Observe that since O commutes with the time derivative, in many relevant cases O defines 
a symmetry for the evolutionary problem if and only if it is a symmetry for the stationary 
one, i. e., if and only if E{(f)) ~ EiOcf) for all states 0, where E is the energy functional. 

In the case of network equations, a new class of symmetries arise in a natural way: the class 
of proportions respected pointwise by physical quantities (e.g., temperature, densities, wave 
functions...) along the edges of a network during the time evolution of a physical process. 
To fix the ideas, consider a closed linear Y subspace of C™ (m being the number of edges in 
the considered network). Then a linear closed subspace of the state space := {L^{Q, 1))™ 
can be naturally constructed as 

y:={feX^: f{x) e Y for a. e. xe (0,1)} . 

We say that V reflects a symmetry of the network diffusion equation if the solution u(-, /) 
to the problem with initial value / satisfies 

vu(tj) = u{t,rf), t>o, 

cf. Definition 15.21 below, where this is formulated in terms of the strongly continuous semi- 
group [e^^)t>o generated by the operator A. In Section [5?^ we will show that this is the case 
if and only if the orthogonal projection V onto 3^ commutes with the operators of the semi- 
group that governs the parabolic problem. We will also show that in the self-adjoint case 
this is equivalent to the fact that £(</>) — C{e^'''^(p) for all s G R, where C is the Lagrangian 
of the Schrodinger system corresponding to the parabolic problem. In other words, we will 
see that V reflects a symmetry for the parabolic problem if and only if it generates a group 
of symmetries for the Schrodinger system. In this sense {e^'^'^) can be considered as an 
equivalent of a gauge group for our dynamical system. We mention that related notions of 
symmetries on quantum graphs have been discussed by several authors, cf. [TT]-|28|-[7]. 

Throughout this paper we will consider directed graphs. This may be disorienting at first, 
since we are always concerned with isotropic physical processes. In fact, all results about 
well-posedness as well as all those concerning positivity and asymptotics of solutions do not 
depend on the chosen orientation of the graph underlying the network, as it can be expected 
(and as it is proved in Section 3). However, we will see in Section 4 that symmetry results 
do in general depend on orientation: in fact, each orientation of the graph corresponds to 
different symmetries. 

We will explicitely consider parabolic systems of diffusion equations in the most part 
of this paper. However, we will discuss in Section 15.21 how symmetry properties of both 
parabolic and Schrodinger problems can be related by means of the theory developed in 
Section 4, see [6, for more details. 



SYMMETRIES IN STRONGLY COUPLED NETWORK EQUATIONS 



3 



2. WELL-POSEDNESS of THE NETWORK EQUATION 

The basic objects we will consider are finite directed graphs, i.e., quadruples G of the 
form (V, E, So, Si) where V = {vi, . . . , v„} and E — {ei, . . . , 6^} are finite disjoint sets and 
So, Si : E — > V are mappings. They associate to an edge e two vertices e(0) := (5o(e) and 
e(l) := Si{e), which are called initial and terminal endpoint ofe, respectively. This promptly 
leads to introducing two matrices X+ = (i^^ ) and T~ = [li^- ) that fully describe the structure 
of the graph. They are defined by 

(2.1) .t, ■■-{]; ^^;^(«).= ^^-' and '{'^^^X^''^^ 

'^■J [ 0, otherwise, '^^ y 0, otherwise. 

Observe that if X := Z+ — is the incidence matrix of the directed graph G as commonly 
considered in graph theory. If \bkj \ = 1, then the edge is said to be incident to the vertex 
Vfc. We define 

r+(vfc) {] e {1, . . . , m} : e,(0) - v^} and r-(vfe) := {j G {1, . . . ,m} : e,(l) = Vfe} , 
and by r(vfc) := r+(vfc) U r^(vfc) we denote the set of indices of all edges that are incident 

to Vfc. 

If there exists an edge e £ E such that either e(0) — and e(l) = v^, or e(0) — and 
e(l) — Vfc, then the vertices v^, Vf are said to be adjacent. Similarly, we say that edges e^, e^- 
are adjacent if there exists a vertex which they are both incident to, i. e., if there exists v G V 
such that 6^(0) = v or 6^(1) — v, and such that 6^(0) = v or ej(l) = v. 

Additionally, we assign to the graph a metric structure that allows us to treat it as a one- 
dimensional manifold and, eventually, to consider partial differential equations describing 
evolution processes taking place on it. Throughout this paper we will always call network 
any directed graph endowed with such a metric structure. A similar if not identical approach, 
based on von Below's theory of C^-networks, has been presented in [29 . 

More precisely, each edge of the graph will be thought of as an interval. For the sake of 
consistency with the notation introduced in (|2.ip . such intervals are parametrized in such a 
way that they have length 1. Whenever we consider a square integrable function / acting 
on the graph G, we may equivalently think of / as a complex- valued function G ^ C defined 
almost everywhere (with respect to the 1-dimensional Lebesgue measure) on the edges of the 
graph, or equivalently as a vector-valued function (0, 1) C™. In this case we will denote 
/ by (/i, . . . , fm)^, where each fj G L^{0, 1) is a function on Cj, i = 1, • . . , rn. Whenever 
point evaluations of / are well-defined, we define with an abuse of notation fj(yk) ■= /j(0) 
if i+ = 1, and /j(vfc) := fj{l) if i^T^. = 1. 

As already emphasized in Section 1, in contrast with the setting which is usual in the 
literature on network evolution equations, we discuss a general model and allow for (possibly 
non-mutual) interactions of non-adjacent pairs of edges, too. The influence of the process 
taking place along the edge onto that taking place along will be descrived by the 
ephaptic coupling coefficient Cji. Such a coefficient is seen as a function on the edge e^: with 
the same convention as above we thus denote Cji(vf) :— Cji(O) or Cji{\ii) :— Cjj(l) if = 1 
or = 1, respectively. 

While the dynamics of our system is described by the coupled diffusion equations in (jl.ip , 
we still have to equip it with suitable conditions in the nodes. To this aim, we introduce 
two tensors defined by 

3+ := Z+ (g) Z+ and 3" ■.= X^ ®X~ . 

We call 3 := 3^ — 3^ the ephaptic incidence tensor of G. Here ® stands for the usual 
Kronecker product of two m x n matrices, defined by {A (g) i?)^f := Ukj ■ bu- We denote by 
entries of 3, 3+, 3_, respectively. In other words, t^/ represents the influence 
of the vertex Vf as an endpoint of on the vertex as an endpoint of . By construction, 
such influences are symmetric, i. e., = i^j for all i, j = 1, . . . , m and all fc, = 1, . . . , n. 
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Solutions of our network diffusion problem have to be continuous in the vertices, i. e., 

(2.2) Ui{vk) = Uj(vfe) for aU i,j e r(vfc), fc = 1, . . . ,n. 

Because of the continuity condition expressed in the equation (|2.2p . we can and will denote 
by the joint value of the components of the vector- valued function u at the node v^. 

Furthermore, we allow (possibly non-adjacent) vertices of the graph to influence each 
other. A natural interaction condition can be formulated as 

m n 

Here, the weighted incidence tensor W :— (w^f), for i, j = 1, . . . , m and l,k — \, . . . ^n, is 
defined by 



UJ 



In fact, in a fashion similar to that considered in [23] we allow for even more general, 
non-local Kirchhoff-type conditions. Such generalized conditions are given by 



m n 



(2.3) ^ ^w^*;fu-(t, v<;) = ^mfe^rf^, fc=l,...,m. 

i,j=l 1=1 1=1 

Summing up, we investigate the strongly coupled system of initial-boundary value diffu- 
sion problems 



(2.4) 



u 



3{t,x) = E(cjiU-)'(t,a;), X e (0, l),t > 0, j = 1, . . . ,TO, 



=1 



Ui{t,Vk) = Uj{t,Vk) dl{t), t > 0, i, j e r(vfc), fc = 1, . . . ,n, 

n in n 



'li S » 

=1 i,j=i e=i 

Uj{0,x) = Ujo{x), X e (0, 1), j = 1, . 



Remark 2.1. By definition of W, whenever C{x) = Id and M = (i.e., if only local, 
synaptic interaction occurs), (|2.3p reduces to the usual Kirchhoff condition prescribing that 
in each node V£ incoming and outgoing heat fluxes agree. 

We introduce := {L^{0, 1))™, which is a Hilbert space with respect to the canonical 
inner product 



™ pi 

{f\9)H^y2 fj{x)gj{x)dx, f,geV. 

:,=1 Jo 





We also consider its dense subspace 

V := {/ e {H\0, 1))" : 3rf^ € C" s. t. {1+^ = /(I), {J-y = /(O)}, 

the space of all if^-functions that are continuous in the nodes of the graph. The subspace 
is a Hilbert space with respect to the canonical inner product 
m „i 

(/ I 9)v ■■= E / {ni^)W^ + Mx)'^) dx, f,geV. 
i=i Jo 

Observe that V is densely and compactly imbedded into X'^, since (C^(0, 1))™ C V C 
(i'(0,l))". 

For the sake of later reference, we recall that a complex (possibly nonsymmetric) matrix 
M = {rriij) is called accretive (resp., dissipative) if there exists /i > such that Re(iV/^|^) > 
(resp., Re(Aff |i^) < — /^|^P) for all e C". We call M positive definite (resp., negative 
definite) if it is accretive (resp., dissipative) and moreover /i can be chosen > 0. 

Throughout the remainder of this paper we will always assume the following. 
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Assumption 2.2. The coefficients Cij are functions of class C^[0, 1]. The matrix C{x) = 
(cy (x)) is positive definite, uniformly on the interval [0, 1], i. e., there exists fi > such that 



Re{C{x)v I v) 



Re ^ c. 



for allx e [0,1], V eC' 



Observe that Assumption 12.21 is weaker than Assum. 2.3]. 

Let us now introduce the KirchhofF operators <i>^ : {H^{Q, 1))™ — > C" defined by 



i,j = l 1=1 



$ U 



V,i=i fci 

and a differential operator on by 



i,j = l 1=1 



\id=i e=i 



(2.5) 




dx ("^l™ da; ) 



— fr — ) 



''dx> 



with domain 

(2.6) D{A) := {/ e (ij2(0, 1))" n V : $+/ - = A/d^} , 

for the matrix M = (rnkh) introduced in (|2.3p . Since I?(A) C V", functions in -D(A) are 
continuous in the nodes. 

With the aim of pursuing a variational approach to our problem, we introduce a densely 
defined sesquilinear form a defined by 

m \ n 

(2.7) a(/,5) {C f \ g')x^ - {Md^ \ ds)c^ = E / c,,{x)f'^{x)'^dx - ^ m^.d/df 

i.j=ll k,l=l 

for f,g^V, which will be later shown to be related to the operator A. 

Theorem 2.3. The operator associated with the form a generates a compact, analytic semi- 
group on . This semigroup is contractive (hence asymptotically almost periodic, too) if 
M is dissipative. If M is dissipative, then the semigroup is strongly stable if and only if 
M*l 0. The semigroup is uniformly exponentially stable if M is negative definite. The 
semigroup is self-adjoint if and only if the matrices C{x), x G [0, 1], and M are self-adjoint. 

Observe that the last result also characterizes well-posedness of the quantum graph asso- 
ciated with ((2^ . 

We stress that if the semigroup is contractive (resp., uniformly exponentially stable), then 
AI is not necessarily dissipative (resp., negative definite), as one sees already in the case of 
a network consisting of a single interval, if one considers the function / defined by f{x) = x 
and M = Id. 

Proof. We show that the sesquilinear form a is continuous and A^-elliptic, i. e., 
• k(/:5)l < ^1 ll/lly WgWv some constant Ki > and all f,g e V, and 



there exist 



> and e R such that Rea(/,/) > - uj\\f\\j^2 for aU / £ V, 



respectively. In fact, the continuity of a is a direct consequence of the Cauchy-Schwarz 
inequality in and of the continuous imbedding of V into (C[0, 1])™, and the constant Ki 
is the maximum over x G [0, 1] of the matrix norm ||C(a;)||. 

In order to prove A^-ellipticity of a, it suffices to observe that (C/' | g')x^ clearly defines 
an X^-elliptic form if (and only if) C{x) is a positive definite matrix for a. e. x G [0,1], 
which is Assumption 12.21 Since there exists K2 > such that 



max \f(x)\ < 

ri;e[0,l] 



/ei?'(0,i). 



6 STEFANO CARDANOBILE, DELIO MUGNOLO, AND ROBIN NITTKA 

cf. [HI Cor. 4.11], it follows that the space of continuous functions over the graph is an interpo- 
lation space between {H^{0, 1))™ and (-^^^(0, 1))™. It then suffices to apply 24, Lemma 2.1] 
in order to treat the lower order perturbation given by {Mdf I fi9)c". Accordingly, by ^ 
Prop. 1.51 and Thm. 1.52] the operator associated with a generates an analytic semigroup 
of angle ^ — arctan Ki . 

Observe that by the Rellich-Kondrachov theorem the embedding of V into is com- 
pact, thus the semigroup is compact. A direct computation shows that a is accretive (i.e., 
Rea(/, /) > for all / S if M is dissipative; and that a is coercive (i. e., it is X^-elliptic 
with w = 0) if M is negative definite. In the first case, the semigroup associated with a is 
contractive, and by [21 Thm. 5.5.6] also asymptotically almost periodic. In the latter case, 
the semigroup is uniformly exponentially stable since the shifted form a — a(-|-)y is accretive. 
Finally, let M be dissipative. Then, by [ini Exa. V.2.23] the semigroup associated with a 
is strongly stable if and only if is not an eigenvalue of the operator associated with the 
adjoint form a*. First of all, observe that if A* f = 0, then necessarily 

f^wf'wh < icr\r)x- = (Md^MOc" < 0, 

thus / is a constant, i.e., / — cl. Observe now that is an eigenvalue of A* (and thus 
necessarily with eigenfunction 1) if and only if 

= {A*l\g) ^ -a* (1,5) = (M*l\dS)c^ for all g e V, 

and since the nodal values of g are arbitrary vectors of C", this is equivalent to saying 
that M*l = 0. Finally, a is self-adjoint if and only if so are the coefficient matrices. □ 



Remark 2.4. It is known that the operator associated with the form a cannot generate 
an analytic, quasicontractive semigroup unless a is X^-elliptic, (see [1] § 5.3.4]), and hence 
unless Assumption 12.21 holds. 



In order to show the well-posedness of our motivating problem, we need to make sure the 
operator associated with a is actually A as introduced in (|2.5p ~ (|2.6p . Having proved this. 
Theorem 12.31 becomes a generation result for A, and in the remainder of this paper we will 
denote by (e*"^)t>o the semigroup introduced above. 

Proposition 2.5. The operator associated with a is {A,D{A)) as defined in (I2.5p - (l2.6p . 

Proof. Denote by {B, D{B)) the operator associated with the form a, which by definition is 
given by 



D{B) := {f eV -.BgeX^ s.t. a{f,h)^{g\h)HyheV}, 
Bf := -g. 



We first show that AcB. Fix f E D{A). Then for allheV 
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Using now the definition of the incidence tensor 3 — J ~ J we can write 



m n 



■m n 



fii^e)'T'ji^k) 

i,j=l k,l=l 



m, n 



n 

e,k=i 



As a consequence 



1 

Thus, for all /i e there exists g — Af e such that 

m \ 

aif,h) = - E / i(^3^f'i)\^)^oi^)<^^ = I 

This completes the proof of the first inclusion. Conversely, let / e D{B). By definition 
there exists g G V such that a{f, h) = — (g | h)^ for all h G V, and accordingly 

m J. n rri \ 

E / Cii(a^)/i(a;)ft.^(a;)dx - ^ miud{d^ = / gi{x)hi{x)dx. 

i,j=l{ e,k=i i=l{ 

Integrating by part the left hand side, we obtain that 

m \ n m n n m \ 

- E / {c,,fj{^)W)dx^TA E E'^^ /^(^^)- E = -E / 3.(2^)^^^^. 

iJ = i Q k=l = l f=l fc,£=l »=1 

which holds for all h E V. In particular, considering /i e (iJo(0, 1))™ vanishing on all but 
one edge of the network, we conclude that 

m 

9i{^) — ''^{cjifjYix) for all x G (0, 1) and alH = 1, . . . , m. 
Similarly, considering h with arbitrary nodal values and arbitrary small X^-norm, we obtain 

m n n 

E fii^f-) ~ E ^'^^'^e =0 for all A: = 1, . . . , n. 

This shows that / G -D(^) and completes the proof. □ 

Having proved analytical well-posedness in an i^-space, one could try to extend this result 
to further L^-spaces, p ^ 2. To this end, a common strategy is to show that the semigroup 
leaves invariant the unit ball of L°° , so that each operator e*'^, t > is contractive on all 
spaces, p G [2, oo], by virtue of Riesz-Thorin interpolation theorem. This has already been 
accomplished in the case of pure synaptic coupling, cf. [HI -[13]. However, we show in the 
following that this approach cannot work in the case of nontrivial ephaptic coupling. 

Theorem 2.6. The following assertions hold. 

(1) The semigroup (e*'^)t>o is real, i. e., it leaves invariant the suhspace of real-valued 
function of X^, if and only if 
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• C{x) e Afm(IR) for all x G [0, 1] and 

• Me Mn{R). 



(2) The semigroup (e*'^)t>o is positive, i. e., it leaves invariant the positive cone of , 
if and only if 

• C{x) is a real valued, diagonal matrix for all x G [0, 1] and 

• the matrix M has real entries that are positive off-diagonal. 

In this case, the semigroup is also irreducible if the graph is connected. 

(3) The semigroup {e*^)t>o is X°° -contractive (resp., -contractive) , i. e., it leaves 
invariant the unit hall of X°" (resp., of X^), if and only if 

• C{x) is a real valued, diagonal matrix for all x G [0, 1] and 

• the matrix M satisfies Re mkk+J2h^k l^fe/il - ^ (resp., Remkk+J2h=ik I'^hk] < 
) for all k = 1, . . . , n. 

Proof. As shown in the proof of Theorem l2.3l the form a is densely defined, continuous, and 
X^-elliptic. Thus, by |261 Prop. 2.5], and by a simple rescaling argument, the semigroup 
(e*")t>o is real if and only if Re / S F and a(Re/, Rg/) G IR for all f & V. Thus, an easy 
computation shows that reality of the coefficients C, M is sufficient. 

Conversely, assume (e*"^)t>o to be real. Let /o G Hq{0,1) real valued and such that its 
support of fj agrees with [a,b] C (0, 1). Define / as a function such that fi = ifo, fj = fo, 
and all further coordinates vanish. By the above characterization of real semigroups one 
has a(Re/,Rg/) — (x)|/o(a;)|^dx G R. Since this construction can be repeated for 

arbitrary a, 6 and we deduce that Cij{x) is a real number for all x G (0,1), and by 
continuity also for all x G [0, 1]. 

Let now f G V such that d{ — 1 and dj. — i. If / vanishes in all further nodes, 
a(Re/, Rg/) = (C(Re/)' | (Rg/)')^^ — mke- As shown above, C(x) is a real matrix for all 
x G [0, 1] and therefore (C(Re/)' | (Rg/)')^^ G R. Thus, niki G R for all fc, ^ = 1, . . . , n. 

In a similar fashion and taking into account [23) Thm. 3.5] and [9] Prop. 3.6], one can 
prove the claimed characterizations of positivity, A'°°-contractivity and, by duality, X^- 
contractivity of (e*'^)t>o. □ 



Additional properties of boundary regularity of solutions of (12. 4p can be deduced by the 
fact that the analytic semigroup operators e*'^ map X^ into [j^^i D{A'^) for all t > 0. 



Proposition 2.7. If u is the solution to (|2.4p . the following assertions hold. 

(1) i( 

Ciju'jY is continuous in the nodes and satisfies a Kirchhoff law, i. 



^(c,,u:)'(t,v,) = ^(c,,u:)'(t,v,) t > 0, e r(v,), £ = 1, 

i,=i 



J2^kidrHt) = I] I]'^''(c,X)"(^,v,), i>0, fc = l,...,n. 

i=l t,i,i=l i=l 

(2) If furthermore the coefficients matrix C is diagonal, then u is of class C°° and its 
derivatives of even and odd order satisfy for all N eH 

uf''\t,y,) = uf^)(t,v,)=:df"'(0, i>0, z,jGr(v,), £=l,...,n. 



m n 



Y^mkedf"\t) = Y.T.^'>r^'\t^^')^ t>0,k=l,...,n. 

3. Symmetry Properties 

In this section we will characterize invariance of different classes of closed linear subspaces 
of the space X^ under the action of {e^^)t>o- The invariance of a closed subspace under the 
action of a semigroup can be characterized as a direct consequence of a result due to E.-M. 
Ouhabaz, see [26l Thm. 2.2]. For the sake of self-containedness we present it in the form 
we will use in the following. Observe that in the view of [9l Cor. 5.2], the invariance results 
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for susbpaces deduced by means of Theorem 13. II can be directly extended to a large class of 
nonlinear, strip- like subsets of . 

Theorem 3.1. Let a:V y.V he a continuous, elliptic sesquilinear form on a Hilbert 
space H, and consider an orthogonal projection P on H. Then KgV is invariant under the 
action of the semigroup (e*'^)f>o associated with a if and only if 

(1) PV <ZV and 

(2) a{f,g) = for all f eRgrr\V,g e KevV n V. 

A relevant class of subspaces of can be constructed as follows: Let F be a subspace 
of C™ and consider 

(3.1) y:^{feX^: f{x) G Y for a. e. x£ (0,1)} . 

We look for criteria for invariance of the subspace y of X^ under the action of the semigroup 
(e*^)t>o- Denoting by K the orthogonal projection of C™ onto Y, the orthogonal projection 
Vk of X-^ onto y satisfies 

(3.2) {VKf) (x) = K (fix)) for a.e. x € (0, 1). 

The aim of this section is to discuss problems that are similar to that presented in the 
following, which also shows an intuitive relation between invariance and symmetry proper- 
ties. 

Example 3.2. Consider a graph G consisting of two edges, both outgoing from a common 
vertex v, i. e., an outbound star. Let C = Id and M = 0. Then the form a is associated with 
the Laplacian with a Kirchhoff condition in Vi and Neumann conditions in the boundary 
nodes. Do initial data that are symmetric with respect to v give rise to solutions to the 
diffusion problem that are also symmetric with respect tov? We can reformulate this question 
and ask whether the closed linear subspace y := {f & X^ : /i = /2} is invariant under the 
action of the semigroup (e*^)f>o. In fact, y = KgVK, where K is the 2x2 matrix whose 
entries equal i. 

Let us reformulate the criterion in Theorem 13. II in our special case. After rewriting the 
form a as a{f,g) = {C f | g')x'^ ~ {Md^ \ d^)cn. observe that the denseness of 14 := {/ G 
V : d^ = x} in X^ for each a; £ C" implies that the condition (2) of Theorem 13.11 holds if 
and only if 

(3.3) (C/' I g')x^ =0 for aU / e Rg -P n V", 5 e Ker V r\V 
and 

(3.4) {Mdf I (is)cn = for aU / e Rg "P n g £ Ker 7? n V. 

We will refer to condition (1) of Theorem l3.1l as to the admissibility of the projection Vk (or 
sometimes of K)^ and to the condition p.3p and (|3.4p as the orthogonality condition with 
respect to Vk of the coefficient matrix C and of the matrix M, respectively. Characterizing 
admissibility and orthogonality is aim of the following subsections. 

3.1. Admissibility. In particular, 'K&cVk and JigVx are isomorphic to (i^(0, 1))*^ and 
(L^(0, l))"", respectively. 

We will now investigate the admissibility of projections of the type Vk in terms of the 
matrix K and of (the incidence matrix I of) the graph G. Let us fix some notation. For 
-4 C {1, . . . , m} we define the vector 

(3.5) 1a ■= (ai)i=i,...,m, where := 
and write t := t_A in the special case of ^ = {1, . . . , to}. 

Lemma 3.3. Let the graph G be connected and the projection Vk be admissible. Then 1 is 
an eigenvector of K. 
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Proof. By hypothesis VkV C V. Consider the function and observe that 

■PkI(x) = Kl and 1 e This shows that on each edge Vk^ is a constant function, and 
since Pk^ G V ail these constants coincide, hence Pk^ = al for an a G C. □ 

Remark 3.4. Observe that Kl e {0,1}, since the only eigenvalues of an orthogonal pro- 
jection are and 1, and that 1 G Ker(Id —K) if 1 e RgK. Moreover, K is admissible if and 
only if Id —K is admissible. Therefore we may assume 1 € Rg K without loss of generality. 



Lemma 13.31 can be used to investigate the invariance of subgraphs. 



Example 3.5. If the graph G is connected, then there exists no proper subgraph G' of G 
such that the linear subspace y {f € : f\G' = 0} of the functions vanishing on G' is 
invariant under the action of {e^^)t>Q- 

Without loss of gererality we may assume that the subgraph G' corresponds to the edges 
e„,/+i, . . . , e„i. The projection onto Y is given by Pk, where 

Of course, 1 is not an eigenvector of K. This result is independent of the matrices C and 
M. 

To characterize admissibility of projections having 1 as an eigenvector we introduce some 
additional notation. We define the 2m x n matrix I and the 2m x 2m matrix K as 

(3.6) I:^{J+,J-V ={fj^Jy^ and K := {^^ 

Observe that K is an orthogonal projection of C^™. 

Lemma 3.6. Let the matrix K he an orthogonal projection of and the let the set Y be a 
linear subspace o/C*. Then the following assertions are equivalent. 

(a) KY C Y; 

(b) y = Kcrxnr ©Rgxnr. 

Proof, "(b) => (a)". Let m e F, i. e., u ^ ui + U2, where ui G KeiKCiY and U2 G RgKCiY. 
Then Ku = Kui + Ku2 = U2 G Y , which proves the claim. 

"(a) => (b)". Let = {b} : i = 1, . . . , tq} be a basis of Keri^ nF and = {&f : 
i = 1, . . . , qo} be a basis of KgK n Y. Extend B^ and B^ to a basis of KeiK and KgK, 
respectively, denoted by 

and 



B'* = B'u{bl ■.i = ro + l,...,r}, 
B^* = B^U{b-':j=qo + l,...,q}. 



J 

Observe that = Ker K®RgK since K is an orthogonal projection. Let u eY . Then 

r q 
i=l 1=1 

with imiquely determined coefficients ai, (3j,i = 1, . . . , r, j = 1, . . . ,q. Now 

ifu~^/3,62= J2 ^3^b]eRgKnY, 

i=l i=qo + l 

since Ku G K by assumption, and hence /Sj = 0,j = go + 1, ■ • • ,9 by definition of _B^*. 
Analoguously it can be shown that ai — 0, i — ro + I, ... ,r hy considering (Id —K)u. This 
shows u G Ker K f]Y ®RgK r\Y. □ 

Proposition 3.7. // the graph G is connected, then the following assertions are equivalent. 

(a) The projection Pk is admissible. 

(b) The range of I is invariant under K, i. e., KRgX C Rgl. 
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(c) There exists a basis o/RgX consisting of eigenvectors of K . 

Proof. We start by proving the equivalence of (jaj) and ^ . Recall that for every f E V there 
exists a vector G C" such that 

The admissibility of the projection is equivalent to the fact that for every f E V there exists 
a vector dPf G C" such that 

[I+Yd^'-^f = VKfiO) = Kf{0), {I-ydP'<f = VKf{l) - Kf{l). 

Inserting the first equation into the second and observing that for all u G C" there exists 
a function / S (iJ^(0, 1))™ which is continuous in the nodes such that d^ — u one obtains 
that (jaj) is equivalent to the fact that for all it 6 C" there exists w S C" such that 

{I+)'^v = K{I+)^u, {I-yv^K{I^)'^u, 

which can equivalently be stated as 

^RgX C RgX. 

The first equivalence is now proved. To see the second equivalence, observe first that the 
existence of the claimed basis is equivalent to RgX being decomposable into RgX = (Kei KCi 
RgZ)©(Rgi^nRgX). Now one can applv the Lemma l3.6] setting Y :— RgXandif := K. □ 

Lemma 3.8. Consider a decomposition G = Gi U G2 into subgraphs such that every node is 
contained either in Gi or G2. On Gi, fix a non-admissible orthogonal projection Vkx- Then 
the projection Vk on G defined by 



K -.^ 

is not admissible. 



Ki 
Id 



Proof. Since Vki is not admissible, there exists a function / G Vi such that Vxif ^ Vi, 
i.e., such that the continuity condition is violated in a node Vfc„. It is possible to extend the 
function / to a function / on the whole graph, such that d-^ = in all nodes of G2. Then 
the function 'Pk/ does not satisfy the continuity condition in v^^, either. □ 

3.2. Orthogonality condition — the matrix C. We are now going to characterize the 
coefficient matrices C which satify the orthogonality condition; in fact, we will show that 
the orthogonality condition is equivalent to the invariance of the range of Vk under the 
coefficient matrix C. 



Proposition 3.9. Let the sesquilinear form a on X be defined as in (|2.7p . with M ~ 0. 
Then the following assertions are equivalent. 

(a) The matrix C satisfies the orthogonality condition p.3j) with respect to Vk- 

(b) The range of K is invariant under the action of C{x) for all x, i. e., 

(3.7) C{x)RgK dRgK for all x 

Proof. Since the space can be decomposed into X'^ = JigVK ©Rg(Id — Pk), the orthog- 
onality condition (|3.3p is equivalent to a{'PK'u, (Id— 'Pi<-)w) = for all u,v V. Using the 
linearity of the derivative and the self-adjointness of the orthogonal projection K, one can 
compute 



a{VKU,{U-VK)v) = [ 

Jo 



{C{x)Ku'{x) I {Id-K)v'{x))dx 
{{ld-K)C{x)Ku'{x) I v'{x))dx, 







where the inner product is the standard inner product in C™. By a localization argument 
{{Id -K)C{x)Ku'{x) I v'{x))dx = Oholdsforevery w,u e y if andonly if (Id-if)C(x)is: = 
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for all X e [0, 1], i. e., C{x)K — KC{x)K for all x £ [0, 1]. Since _ftr is a projection this is 
equivalent to condition (|3.7p . □ 

3.3. Orthogonality condition — the matrix M. Next we characterize the orthogonality 
condition for the matrix Af , i. e., we want to find equivalent conditions to p.4p . where (|3.4p 
can alternatively be stated as 

(3.8) (Md^^^ I -T^'<)a) ^ for aU f.geV. 

If it is satisfied, we will say that the matrix M satisfies the orthogonality condition with 
respect to K, frequently omitting any reference to K . For these investigations we introduce 
the matrix 

where we denote D the diagonal matrix with diagonal entries |r(vfc)|, the degrees of the 
nodes v^. Please note that the matrix Ai only depends on AI and on the graph structure, 
but does not depend on the orthogonal projection K. 

Lemma 3.10. If the graph G has no isolated nodes, then the following assertions hold. 

(1) df = i?-ijT(/(0),/(l))T for every f€V. 

(2) RgZ={(/(0),/(l))TeC2-:/GT/}. 

Proof. First we will prove the formula 

(3.9) X+ {I+Y = diag(r+(vfe))fc=i,...,„. 
In fact. 



Since each edge originates from exactly one node, we obtain that — for all k ^ I. 

Thus, 



0, otherwise. 



1=1 

Since T'^- equals 1 exactly r+(vi.) times and equals otherwise, the proof of formula p.9p is 
complete. The analogous formula (X^)^ ~ diag(r^ (vfe))fe^i_..._„ can be proved likewise. 
As a consequence, we obtain 

D = i+ {i+y +1- {j-y = i^i. 

To prove (P), let / G By definition, there exists d^ € C" such that 
(3.10) id^ 



ffiO) 



We show that the vector D \^f \\)j satisfies the condition l|3.10p as well. A direct 
computation shows that 

iD-^r =iD-^ndf =iD-'Ddf =idf = (^j^jjj^ . 

By the uniqueness of d^ , the proof is complete. 
For dH) notice that since {df : f eV] ^ , 

Jigi^[iv : V e c"} = {idf : / e y} = --f^v 

This completes the proof. □ 
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Proposition 3.11. Assume the orthogonal projection Vk to be admissible. Then the matrix 
M satisfies the orthogonality condition (|3.8p if and only if 



(3.11) RgMKI cRgK. 

Proof. We will use the orthogonality condition as stated in (|3.8p . By Lemma [3.101 (1). one 
obtains for all f,g € V 

By Lemma |3. 101 (2). the orthogonality condition is then equivalent to 

(3.12) {MD-^J^kv I D-^i'^ {Id -k)w) = 0, Vu, w e Rgl. 

Since X has real entries and D^^ and (Id— i^T) are self-adjoint, p.l2p is equivalent to 

((Id ~k)JD-^MD-^i^kv I w) = 0, Vu, w e Rgi, 
i. e., using the definition of M 

((Id -k)Mkv I w) 0, Vw, w e Rgi. 
This can be equivalcntly expressed as 

((Id ~k)MkP^^^v I P^^^w) = 0, V«, w € c™, 
where ^j^gj is the orthogonal projection onto the range of I. Since Pji^j self-adjoint, 
{{M~k)MkP^^iV I P^^^v) = {Pn,xi^d-k)MkP^^iV \ w) 

for all e C"\ In fact, we have just proved that the orthogonality condition (|3.8p is 
equivalent to 

Pnsi{ld^k)MkP^^i^O, 

which is, finally, the same as 

(3.13) PR,iMkP^^i = PR,iKMkP^^i. 

The equation p.l3p is the key to prove the claim. Because of the admissibility of K, 
RgkX C Rgi by Proposition [321 Moreover, one sees by the definition of A4 that Rg C 
RgX, which implies Rg KM C Rgi. Considering both inclusions, one obtains that (|3.13p 
is equivalent to MkP^^j- = kM.kP^^j-. In fact, we are asking that k acts as the identity 

matrix on RgMkP^^j- = RgMki, i. e., the orthogonality condition (|3.8p is equivalent to 
Rg Mki C Rg k. This concludes the proof. □ 

Although it is easy to check those range inclusions numerically for concrete examples, the 
following sufficient conditions may be more convenient in some cases. 

Corollary 3.12. Consider an admissible projection Vk- IfRgMi C Rgk or RgMk C 
Rg-R', then M satifies the orthogonality condition. 

Proof. To see that RgA^Z C Rgk is sufficient, observe that admissibility of Vk implies 
RgkJ C Rgi, and hence condition (|3.1ip is fulfilled. Moreover, RgMK C RgK is also 
sufficient by a similar argument, since RgKX C RgK. □ 

As an application we give a a simpler characterization for a special subspace of C" in the 
case of a bipartite graph. More precisely, we consider the smallest subspace of C™ whose 
orthogonal projection K satisfies Kl — 1, i. e., Y — {(c, c, . . . , c)'^ | c e C}. 

Proposition 3.13. Let the graph G be bipartite. Then M satisfies the orthogonality condi- 
tion with respect to the orthogonal projection K := (^)^ ^ if and only if there exist 
values (ay )ij=i.2 such that 

aii|r(v£)| = X^feii "^^fc: (^i2\T(ye)\ =Ylk=m+i'^'ik for all £ = I, . . . ,ni, and 
Q!2i|r(v<;)| = X^fcii Wtt, a22|r(vf)| = X]fe=„i+i for all£^ni + l,...,n. 
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Proof. First we show Rg KX = Rg K, and then we prove a characterization of those matrices 
Ai that leave invariant the subspace KgK. Then the claim will follow by Proposition [STTlJ 
Without loss of generality we may assume that for k = 1, . . . ,ni the node Vj; has only 
outgoing edges, i. e., r(vfe) = r+(vi,), and that for A: = ni + 1, . . . , n the node Vfc has only 
incoming edges, i. e., r(vfe) — (vfe), reordering the nodes otherwise. We are going to prove 

(3.14) Rg/fX-RgX- 

The second equality follows from the definition of K. Moreover, RgKI C KgK is obvious. 
Since there exists f € V such that — Lemma [3. 101 (2) implies G R-g^^- 

Analoguously = l[ni+i,...,n} yields l{m+i,...,2m} G R-g^- We have already observed that 
l{i,...,m} s-iid l{m+i,...,m} ^^c invariant under K, which implies 

(l{l,...,m}, l{m+l,...,m}) C RgKI, 

and this proves the claim p.l4p . 

Next we characterize the matrices M which leave Rgi^ invariant. For this we use the 
bipartite decomposition of 

M'" := D-^MD-^ = 



|rK)||r(v,)|,,,^,^ 

which is induced by the bipartite decomposition of the graph G, i. e., we write 



Ar 



where Mf^ e M„,^„,, M^^ e M„,,„_„,, M:^^ € M„_„,,„,, and M^^ € M„_„,,„_„,. 
Moreover, since G is bipartite, the incidence matrices decompose into 

andX-^^O 



; v^2" 

where X+ € M„^,m and I2 e Mn^m.m- 

We will use the above decompositions in order to obtain a useful formula for Ai. Using 
the definition of Z, as formulated in (13. 6|) , we first compute 

\I- M'"I+ I- M'^I- j 

which yields, inserting the decompositions, 

M 



^2 ^"^21^1 ^2 ^^22^2 



It is possible to identify the block-matrices appearing in the above expression. In fact, the 
following identities hold. 



It M^^I+ 



-^2 ^^22-^2 



1^ Mi2^2 

ifM^^lt 



™ei (0)63(0) 



|r(e,(0)||r(e,(0))|,^^^.^,^ 

"le.(l)e,(l) 



|r(e,(l)||r(e,(l))| 

™ei (0)63(1) 



|r(e,(0)||r(e,(l))|,^^^.^^^ 

™6i(l)6j(0) 



|r(e,(l)||r(e,.(0))| 



Here we write TOe;(o)ej(o) '■— ^ki if ei(0) — v/j, 6^(0) = v^, and analogously for ej(l) 
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We have already observed that {l{i,...,m}, l{m+i,...,2m}} is a basis of RgK. Using the 
above decompositions, one can compute 



(l+ A/{"iI+^lc"A _ f{^j=i |r(e.(o))||r(e,(o))i; .^^^^^^ ,^ 

[I2 A/2"iXi+) Ic'-y |r(e.(llli[ne"'( 



I 1 I n I "J / I I II I I 1 \ - -I , 

.(o))l^=i,...,. 

In these sums each ej(0) appears exactly |r+(ej(0))| times. Collecting the same summands 
we may write 



A^l{l....,m} 



2^k=i\^ i^^fe^l |r(ei(i))||rK) 



In fact, since |r+(vfc)| appears as a factor, only for the first ni vertices the summand does 
not vanish. Thus, we see that 




(Efeiilr (vfe)||r(e"o»||r(v,)|),^i^ 



|r(e,(i))||rK)| 

Since |r+(vfe)| = |r(vfe)| for fc = 1, . . . , m, 

Eni "T.e^(o)fc 
k=l |r(e,(0))| 

/ Y~>"1 ™ei(l)fc 

A^'==i Ir(e.(i))l 

We can easily check whether G R-g^- Using (|3.14p one sees that this is the 

case if and only if sums above do not depend on e^, i.e., A4t[i^,,,^rn} S RgA" if and only 
if there exist an, 0:21 G C such that aii|r(v^)| = X^feli "^^fc ^^r all i — l,...,ni, and 
a2i|r(v£)| — X^feli '^ffc f'^'' all = ni + 1, . . .. By a similar computation, one can also see 
that A^l{,n+i,...,2m} G Rg-^ if and only if there exist 0112,022 G C such that Q;i2|r(v£)| = 
Sfc=„i+i "^ek for all £ = 1, ... , ni, and a22|r(v£)| = J2k=ni+i "^^k for all ^ = rii + 1, . . . , n. 
This completes the proof. □ 

Example 3.14. Consider a regular, bipartite graph G, with the bipartite node decompo- 
sition G = {vi, . . . , v„j^} U {v„j+i, . . . , v„}. Set 7i2 '■— n — ni and consider row-stochastic 
matrices Mij G Mn-^nj- Then for arbitrary aij G C all matrices of the form 

^aiiA/11 Q!i2Mi2^ 
^a2lM21 Q!22M22^ 

satisfy the orthogonal condition with respect to K defined as in Proposition 13. 131 



M 



4. Classes of graphs 

In this section we will discuss some classes of graphs, combining the results of the pre- 
ceeding sections. We present some (non-standard) graph theoretical definitions we will use 
through this section. 

Definition 4.1. Let G a graph with no isolated nodes, i. e., such that r(vfe) > 1 for all 
k = 1, . . . ,n. 

• We call the graph G completely unconnected if G is the union of disjoint compact 
intervals, i. e., if G is a regular graph of degree 1. 

• We call the graph G an inbound ( respectively, outbound ) star, if there exists a node 
v/j such that 6^(1) — Vk, (respectively, if ej{0) — Vk), for all j = 1, . . . ,m. We call 
the graph G a star if it is an inbound or outbound star and Vk the center of the star. 

• We call the graph G bipartite if each node has only either incoming or outgoing 
edges. 
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• We call the graph G Eulerian if all nodes have the same number of incoming and 
outgoing edges. 

• We call a graph G a layer graph if there exist disjoint sets Vi, . . . ,Vl such that 

— Qj{0) £ Vp implies ej(l) £ V^+i for all p — 1, . . . , L — 1, and 

— ej(0) G Vl implies 6^(1) € V\. 

Nodes belonging to Vp are said to lie in the p"' layer. Edges outgoing from nodes in 
the p^^ layer are also said to lie in the p*^ layer. 

• We call a layer graph symmetric if the incoming and outgoing degrees of the nodes 
only depends on the layer, i.e., if there exist numbers l{p),0{p) G No such that 
|r+(v)| = I{p), |r^(v)| = 0{p) for all nodes v in the p*^ layer. 

4.1. Bipartite and Euler Graphs. It is possible to characterize some classes of graphs 
by the admissibility of the matrix from Proposition 13.131 

Theorem 4.2. Consider the orthogonal projection K defined by 

(4.1) K 

Then Vk is admissible if and only if G is bipartite or Eulerian. 

Proof. Fix f ^ V and observe that T'a'/ always lies in {H^{0, 1))™ since every component 
is a linear combination of functions. So V is invariant if and only if Vkf is continuous in 
the nodes. Let Vi C V denote the set of all vertices having outgoing edges, and let V2 C V 
denote the set of all vertices having incoming edges. We distiguish two cases. First, assume 
Vl n V2 = 0. Then G is a bipartite graph. 

On the other hand, if Vi n V2 / 0, then by definition oi K a. vector d^''^^ exists if and 
only if 

(4.2) ylM^f^m^ 

We show now that the equality (|4.2p is equivalent to the graph being Eulerian. First, 
assume that (14. 2p holds for every f eV . Fix an arbitrary \Jk €\l and choose f G V such 
that d^ = Iji}. Then 



-|r+K)| = E^ = E^ = -rK)|. 

m ' ' ^ — ' m ^ — ' TO TO ' ' 

Thus it is necessary that |r^(vfc)| — |r+(vfe)| holds for every k ~ l,...,n. Conversely, 
assume that |r~(vfe)| = |r+(vfc)| holds for every fc = 1, . . . , n. Then 

E ^ - - E |r+K.)| 4 = - E |r-K)| 4 = E 

m TO ■'^ ' 'to ■'^ ' ' ■'^ TO 

j=l k=l k=l j=i 

Hence (|4.2I) is satisfied, so this condition is also sufficient. 

It only remains to show that indeed for every bipartite graph K is admissible. To see this, 
note that for an arbitrary f G V the vector d^'*-'-^ can be chosen to equal X^I^Li 
components belonging to nodes in Vi and to equal X^IIli ^'m^^ components belonging 

to V2. This shows continuity of VKf in the nodes, thus implying Vxf G V. □ 

Remark 4.3. The matrix K defined in (|4.ip acts on a vector v G C™ by substituting each 
component by the average of all components of the vector. The range of such a matrix is 
thus one-dimensional, and one sees that 

RgVK ^{feV:f, = f, for all f, j = 1, . . . , to}. 

Such functions are symmetric in the sense that they are equal on each edge at the same point 
of the parametrization. In fact, Theorem 14.21 characterizes the admissibility of projections 
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whose ranges consist of the functions that are symmetric on the network. It thus gives a 
first answer to the problem stated in Remark [ 



4.2. Stars. Main result of this subsection is a characterization of stars in the class of the 
simple graphs. We first investigate the admissibility of projections. 

Proposition 4.4. The following assertions hold. 

(1) The graph G is completely unconnected if and only if Vk is admissihile for all or- 
thogonal projections K . 

(2) Let G be a simple, connected graph. Then G is a star if and only ifVx is admissihile 
for all orthogonal projections K with eigenvector 1. 

Proof. (1) Since the graph G is completely unconnected, the continuity condition in V is 
empty, and therefore each Vk is admissible. Conversely, if G is not completely unconnected, 
then it is possible to decompose G into the disjoint union of a connected graph Gi with 
TOi edges and the remaining graph G2. Let Ki be an orthogonal projection of C™^, which 
does not have 1 as an eigenvector. Lemma [3.31 and Lemma [3.81 assert that the orthogonal 
projection 

(Ki 0' 
V Id, 

is not admissible. 

(2) Without loss of generality, we prove the claim for an outgoing star with center Vi 
and with the natural numbering of the other nodes. Let the graph G be a star and K he a, 
projection such that Kl = 1. In fact, for this star 



1 = 



1 
Id„ 



Since now K has 1 as eigenvector to the eigenvalue 1, one can compute 

f Kl Q \ ft 0' 



KI — , , — , 

V KldmJ \0 K 

It is now clear that Rgi^T C RgX, and this implies the admissibility of Vk. Conversely, 
assume that the graph G is not a star. One sees that this implies the existence of an 
undirected path of length 3. We will denote it by 61,62,63, possibly relabeUing the edges. 
Our strategy is the following: for each graph that is a path consisting of 3 edges we construct 
a non-admissible projection Vl where LI = 1. We then consider the projection Vk, where 
K is 

Then, by Lemma [3.8i we conclude that Vk is not admissible, although 1 s an eigenvector 
oiK. 

First, consider cycles of length 3. Since each edge can be directed arbitrarily, there are 8 
such graphs. Let us start with the case of a not strongly connected graph. Such graphs are 
neither Eulerian nor bipartite. Thus, Theorem 14 . 21 provides an example of an L as requested. 
If the graph is a (directed) cycle such that ei(0) = vi, consider the projection 




and the function / defined by f{x) (x, 1 — a;, 0)^ £ V. One sees that f eV but "Pa'Z ^ V, 
since VKfix) = (5, i, 0)^ for a. e. x £ (0, 1). 

Consider now the lines of length 3. We split this into three possible cases: G may be 
bipartite line, a (directed) line, or neither a (directed) line nor a bipartite graph. In the last 
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two cases the graphs is neither bipartite nor Eulerian, and hence we can use Theorem 
again. In the case of a bipartite hne, let us consider the projection 






for the parametrization ei(0) — Vi, ei(l) — 62(1) = V2, 62(0) = 63(0) — V3, and 63(1) — V4. 
Consider the function j{x) := (a;,a;,0)^. Again, j but Vk! ^ since VKf{x) — 
(|, j;, |)^ for a. e. x G (0, 1). This completes the proof. □ 

Remark 4.5. In Proposition 14.41 (2) we have assumed the graph G to have no multiple 
edges. In fact, it is not possible to relax this condition, since all orthogonal projections with 
eigenvector 1 are admissible on all connected graphs consisting of 2 nodes and m edges for 
each m G N and each orientation of the edges. 

Now we investigate the orthogonality condition for diagonal matrices C. This will show 
that for a wide class of matrices C there cannot exist non-trivial invariant subspaces of the 
form considered in this paper. 

Lemma 4.6. Let D be a constant diagonal matrix with entries di > 0. Then the following 
assertions hold. 

(1) Assume the coefficients di to he pairwise different. If K is an orthogonal projection 
with eigenvector 1 such that VRgK C Rgi^, then K is trivial, i. e., K — Id or 
K = 0. 

(2) Assume that there exists ia ^ Jq such that di^ =^ dj^. Then there exists a nontrivial 
orthogonal projection K with eigenvector 1 such that PRgX C Rgif . 

Proof. Observe that it is possible to compute the powers of T> explicitly, since it is diagonal. 
In fact, V'^ = diag(cf for every k e Nq. 

(1) Since K is an orthogonal projection and 1 is an eigenvector, either Kl = 1 or Kt = 0. 
If Kl = 1, i. e., 1 £ Kg K , we see by induction that {d'l^d^, . . . ,dt) = T^'^t eRgK for every 
/c G N since KgK is invariant under the action of T>. Now, the matrix V := ((iij)ij=i,...,m, 
defined by 



dij := dl i,j = 1, 



, m 



is the Vandermonde matrix induced by the vector ((ii)i=i....,m, which is regular since the di 
are pairwise different. From this we see RgK = C'^, i. e., K ~ Id. 

If on the other hand Kl — 0, then fix w G RgK. Since 1 is in the kernel of K, 
Rg-ftT C (1)^. Since the range of K is invariant under the action of the matrix C, we 
obtain (C"u | 1) = di'Vi = for every n G N. In particular, v satisfies the equation 

Since V is regular, we obtain v = 0, which implies Rgif = {0}, hence K — 0. 

(2) In order to prove the second assertion, let i j such that di — dj. Consider 

Y := span{P"l = «, . . . , d;;) | n G N} C C", 

and let K be the orthogonal projection onto Y. Since 1 G F, Kl = 1, and in particular 
K 0. Moreover, Vi = vj for all u G F, which implies KgK ^ C™. As a consequence 
K ^ Id. Finally, T>Y C Y, and hence the range of K is invariant under the action of V. 
This completes the proof. □ 

Combining the previous two statements we deduce the following. 

Proposition 4.7. Let the graph G be connected. If the coefficient matrix C is diagonal and 
M = 0, then the following assertions hold. 

(1) Let C he constant and Q he a star. If there exist iq , jo such that Cig — Cj^ , then the 
suhspace 

y ■■= {f e I f,„{x) = fj„ix) for a. e. x G (0, 1)} 
is invariant under the action of (e*'^)t>o. 



SYMMETRIES IN STRONGLY COUPLED NETWORK EQUATIONS 



19 



(2) Assume the coefficients Ci{xQ) to be pairwise different for some Xq G [0,1]. If Y is 
a nontrivial linear subspace of , then y, defined as above, is not invariant under 
the action o/(e*"^)t>o- 

Proof. (1) Without loss of generality, assume ci = C2. Consider the subspace 

Y := {v e : Vi = V2} 

and let K be the orthogonal projection onto Y. Since 1 e F, Kl = 1. Furthermore, 
by Proposition 14.41 K is admissible, since G is a star. Let v E Y. Computing now Cv = 
('^i^j)i=i,- -,ni shows that Cv € Y, since ci — C2 and vi = V2. This shows that CKgK C 
Rg K, thus completing the proof of the first claim. 

(2) Let Y C C™ be a linear subspace, and let K be the orthogonal projection onto Y. 
Remember that the invariance of this subspace is equivalent to the fact that K is admissible 
and that the sesquilinear form a satisfies the orthogonality condition with respect to Vk- 
If K is not admissible, then the proof is complete. Thus, assume that K is admissible. 
Since the graph G is assumed to be connected, 1 is an eigenvector of K. The orthogonality 
condition is equivalent to C{x)Y C Y for all x e [0, 1], according to Proposition 13.91 Using 
Lemma 14.61 we see that since the diagonal entries of C{xo) are pairwise different, this is not 
possible for non-trivial K. Hence the proof is complete. □ 

4.3. Layer Graphs. In this section we prove an admissibility result for symmetric layer 
graphs. We start fixing a canonical numbering of the edges of a layer graph. First observe 
that the node decomposition induces an edge decomposition E — U^^iEp by setting 

E'p := {e G E : e lies in the p"^ layer}. 

After relabeling the edges we may assume that there exist Lp, p — 1, . . . , L + 1 satisfying 

(1) Li =0; 

(2) 8^(0) = ej(0) or 6^(1) = ej(l) implies ip-i < i,j < Lp for some p; 

(3) ei(0) — ej{l) implies ip_i < j < Lp < i < ip+i for some p. 

The numbering obtained in such a way has the property that is in the p**^ layer if and 
only if Lp < « < Lp+i. In fact, all edges such that i < ip+i are in any of the first p layers. 

We are going to exhibit a class of admissible projections. Altough the result is not a 
complete characterization, it is optimal in a sense we will explain later. 

Proposition 4.8. Consider a symmetric layer graph G and the orthogonal projection K 
(4.3) 

where \Ep\, p = \, . . . ,L denotes the number of edges in thep^^ layer. Then Vk is admissible. 

Proof. One has to check the continuity condition for each p — 1,...,L — 1 in every node of 
the p**^ layer. Define the auxiliary function 

A : fc t-* layer of the node v^. 

We thus have to check continuity in those nodes such that X{k) — p, p ~ 1, . . . , L — 1. 
The set X^^{p) can be represented in the form 

X^^ip) = {k:3ie {Lp + l,...,Lp+i} s.t. e,(l) = Vfc}, 

as well as in the form 

A-i(p) = {k:3ie {Lp+i + 1,. ..,Lp+2} s.t. e,(0) = Vfc}, 

whenever the expression is defined. By the definition of K, one sees that for allp — 1, . . . , L — 
1 and all i, j = Lp + 1, . . . , Lp+i the identities 

(4.4) VKf^{l) = VKfjil), rKf^{0) = VKfj{0). 
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hold. As a consequence, for layers having incoming or outgoing degree 0, the continuity is 
obvious. Assume now that I{p) ^ and 0{p) ^ 0. 
For an edge in the p^^ layer and for / £ V, 

i=Lp + l I PI feGA-i(p) ^ 

Recall that since our graph is symmetric, the incidence degree |r~(vfe)| only depends on the 
layer, and therefore we can write 

With analogous computations we obtain for edges in the p + 1 layer 

\0{p)\ 



■p+n 



/ceA-i(p) 



Observe that the identities l-Bp+i] = |A"i(p)||0(p)| and \Ep\ = |A-i(p)||/(p)| imply \0{p)\\E. 



\I{p)\\Ep\ ^. We have thus proved that {VKf)i{^) = CPRDji^) for all 6^,6^- such that 
i G A~^(p), j e X~^{p+ 1). This completes the proof. □ 

Corollary 4.9. Let G be a symmetric layer graph. If M = and C = c(x)Id /or some 
function < c G C^[0, 1], then the space 

y:^{feL^:f,^ f, for all i,j G r\p), p = 1, . . . , L} 

is invariant under the action of {e*'^)t>o- 

Remarks 4.10. 

(1) The class of the layer graphs is not a common object in the graph theoretical literature. 
In fact, layer graphs are nothing but (directed) p-partite graphs, for which collapsing 
the components of the graphs to a single vertex leads to a finite line or to a cycle. In 
particular, homogeneous trees of finite depth are symmetric layer graphs. Such graphs 
play a role in the investigation of biological neural networks. 

(2) The symmetry condition in Proposition 14.81 cannot be relaxed. To see this, consider 
the following simple example. Let G be an outgoing star of order two and consider two 
copies of G. Identifying two of the external nodes defines a layer graph. One can show 
that the orthogonal projection defined in (14. 3p is not admissble, due to the two free 
nodes in the second layer. 

(3) It seems to be possible to extend the result of Proposition 14.81 to non-symmetric layer 
graphs, requiring some weaker condition and suitably weighting the projection of (|4.3p 
according to the degrees. However, such results are quite techical. Presenting them in 
detail goes beyond the scope of this paper. 



p^ 



5. Applications 

5.1. Ephaptic coupling of biological fibers. In the modern neurobiology's early years 
it was common sense that neuron should communicate with each other remotely, only by 
means of their electrical activity. In this context, the theory of so-called ephaptic connection 
was forged in the 1940s by A. Arvanitaki, Nobel laureate B. Katz, and H. O. Schmitt, 
cf. [21 [TB] . Such a theory was thought to be surpassed after the newly invented electron 
microscopes allowed in 1954 to finally prove the existence of chemical synapses. 

Although synaptical connections are ultimately stronger and more common, more recent 
experiments have however found evidence of ephaptic effects in several animals and even in 
human patients. While experiments have been conducted in real neuronal networks, to the 
best of our knowledge mathematical models of ephaptic connections have only been treated 
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in [T3j. Though, there is some hterature for the special case of bundles of (synaptically) 
unconnected nerve fibers of infinite length, cf. |15[ SI [5] and references therein. 

Although the derivation in the quoted articles is different, the mathematical models pre- 
sented in [15l §4], m §4], and [5] are comparable. Possibly up to linearization, they describe 
ephaptic interaction within a myelinated nerve fiber of m axons of infinite lengths by a 
system of diffusion equations of the form 

E;ii(ci,";-(i,-))'(a^), t>0, xeR, 

Er=i(c™j"j(^'-))'(a^), t>0,xeR, 

where Ui{t, x) is the electric potential of the i'^ axon at space x and time t. The coefficients 
{cij) are positive constants that represent the ephaptic effect on the i*'' axon due to the 
activity of the j**^ one. We emphasize that the mutual interactions and therefore the matrix 
{cij) are in general non-symmetric. 

Whenever potential transmission in neuronal networks is mathematically modelled, neu- 
robiologists usually assume that some form of Kirchhoff law holds in the nodes, as well as 
continuity of potential. In the easiest linear case, this amounts to saying that in each node 
the total incoming electric flow equals the total outgoing one, possibly up to some form of 
dissipation, cf. [22]. As we have seen in Remark 1 2. 1[ the natural generalization of Kirchhoff 
node conditions to the case of strongly coupled network equations is given by (|2.3p . 

This motivates us to consider ()2.4|) as a model for transmission of potential in (passive) 
nerve fibers where ephaptic effects hold. The following results allows to easily discuss also 
the computationally hard case of numerous contiguous neurons. 

Proposition 5.1. If the coefficients Cij satsify 

(5.2) c,,>^ '"-^+"^'' , z = l,...,m, 

then the initial value problem associated with (jS.ip is well-posed. 

Proof. By the results of Section 2, the initial value problem is well-posed if the coefficient 
matrix C is coercive. By Gershgorin's Circle Theorem, we directly obtain that (|5.2p implies 
coercivity of the matrix C, and the assertion follows by CoroUarv 12.31 □ 

The coefficients (cij) are phenomenological constants that have to be determined exper- 
imentally. As already observed in [9,, § 4.1], the model proposed in [TS] (i.e., Cy = c for 
all seems to be ill-posed in the light of Remark 12.41 whereas in the models proposed 
in [l]-[5] the possibility to apply CoroUarv 12 . 31 depends on the values given to the coupling 
parameters. 

In all models of ephaptic coupling considered above, the coefficients are assumed to satisfy 
Si!li — consti for all j and X^j^i '^v ~ const2 for all i. Then by Theorem 14.21 one can 
say that a necessary condition for the subspace of pointwise equal functions to be invariant 
under the action of (e*'^)t>o is that the neuronal network is either bipartite or Eulerian. 
In fact, assuming for the sake of simplicity that no dissipation happens in the nodes (i.e., 
M = 0), one deduces that there exists two function Ci,C2 ■ [0,1] C such that for all 
^ G [0; 1] SjLi (a;) = Ci{x) for all i and Cij{x) = C2{x) for all j. 

Observe that by Proposition l2.6l even if the system is governed by a contractive semigroup 
in X'^ (which is the case if M is dissipative), no contractivity property holds with respect 
to the norms || • ||i and || • ||oo unless C is diagonal. In other words, the system's potential 
may increase both globally and locally, as soon as ephaptic effects are actually considered. 

5.2. Quantum graphs. Consider a finite network of thin waveguides ei, . . . , e^, of (possibly 
different) lengths . . ,£m- Discussing the propagation of wave functions, i.e., studying 



(5.1) 



ui{t, x) 

ilrnit, x) 
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the evolution of a system of Schrodinger equations 
dv ■ d'^v ' 

over such hnear structures - usually called quantum graphs - has become a relevant topic 
in recent years, see e.g. and references therein. Kirchhoff or more general 

self-adjoint conditions are usually imposed in the nodes of quantum graphs. 

In order to define an Hamiltonian associated with the quantum graph, observe that after 
a change of coordinates the above equation reads 

= .E(0,l),teR. 

The Hamiltonian is thus given by lA, where (A, D{A)) is the operator introduced in (|2.5p - 
(|2.6p and associated with the form a. Here we are considering coefficients 

I otherwise. 

Thus, the operator A is self-adjoint if and only if the ephaptic coupling and nodal coefficient 
matrices C{x), x £ [0, 1], and M are both sclf-adjoint, which we assume throughout. Then 
by Stones's theorem lA generates a unitary group that governs the evolution on the quantum 
graph. As in classical field theory, we introduce the action functional S for the time evolution 
of the quantum graph (resp., of the parabolic problem), which is defined as 




(resp. 



ipi H 



i.e., S{->p) = /(f (^i{'ip\'ijj)x^ +a{^p,'4})j dt (resp., 5('0) = /J (^{ip\ip)x^ +a{ip,->jj)j dt). Here 

we have imphcitely assumed that ?A £ C^{[0, T]), X^) n C([0, T], V) for an arbitrary T > 0. 
Our aim is to discuss symmetry property of the system, in the following sense. 

Definition 5.2. 

(1) We call a Co-group {U{s))seR onX"^ a symmetry group for the system of Schrodinger 
equations (parabolic equations) over the network if S^ijj) = S{U{s)'iIj) for all s £ R, 
where {U{s)^p){t) := U{s)i:{t), t G [0,T]. 

(2) We say that a hounded linear operator V on reflects a symmetry of the parabolic 
network problem ifVKC^^ — c^^Vk for all t>Q, i. e., if projecting the initial value 
and then studying the corresponding time evolution is equivalent to projecting the 
solution curve of the original problem. 

Since U{s) does not act on the time variable, one sees that due to the time continuity 
oi ip s. self-adjoint bounded linear operator on X'^, i.e., an observable V of the physical 
system, is the infinitesimal generator of a symmetry group (e*''^)^^^ if and only if it satisfies 
a{ip, V) = a{e'''^ip, e*^^V) for all s G R and all V' G V^. 

We consider the case of a closed linear subspace y constructed as in (13. 1|) and discuss 
observables V — Vk given by orthogonal projections of the state space X"^ :— (i^(0, 1))™ 
onto y that satisfy p.2p . A justification for the use of the term "symmetry" in Definition l5.2l 
is given in the following. 

Proposition 5.3. Let Vk be an observable of the system as defined above, and assume K 
to be admissible. The following assertions are equivalent. 

(a) The projection Vk reflects a symmetry of the network parabolic problem. 

(b) The subspace y — Rg'Pif is invariant under (e*'^)t>o. 
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(c) a{VKilJ,^/j) = a{rK'(p,'PK^) for all G V. 

(d) The projection Vk generates a symmetry group of the parabolic network equation, i. e., 
a{ip,ip) = a(e*''P^V,e*^'^'^V') for all s e R and all ip e V . 

(e) The projection Vk generates a symmetry group of the network Schrddinger equation. 

(f) The subspace y — KgVK is invariant under (e'*'^)tgiR. 

Proof. Note that the invariance of 3^ = KgVK under the action of (e*'^)(>o is equivalent to 
(0) VKe'^^VK = e'^^Vx for all t > 0. 

"(jll 0" This is obvious, since Vj^ = Vk- 
"dHl) ^ tlj)" Since Vr and e*^ are self-adjoint, 

Vkc'^ = {c'^VkT = {'PKe'^'VKy = VKe'^VK = e'^Vx- 
"(0 ^ jg)" By Theorem O © is equivalent to a{VKf, (Id -VK)f) = for every 
/ e But this is (jsj). 

"(jSj) ^ ©" Since Vk is a projection, 

OO j OO j 

^.v^ = ^ ^T'^ = ^ + Id = (e^ - \)Vk + Id . 

Using this representation we see that 

+2 Re a ((e'" - \)Vk^. V') + V') 
= a (T'ktA, T'k^A) - 2 Re e'^a (Pi^^A, Vk^^) + a (T'/f^, PkV) 

+2 Re(e*" - l)a (P^^, V') + alV-, V-) 
= 2a{VKi^,VKi^)-2a{VKi^,^)+aiij,i^). 

Thus jdj is equivalent to 2a {Vk^P, Vk^^) — 2a {Vk^P, -0) = for every tp which is |cj). 

"(jdj <^ ([ej)" Both statements are equivalent to a{ip,%p) = a(e**^^^/^, 6*"'^^^/') for aU s G IR 
and all ip & V, since e'*'^^ is an unitary operator that commutes with the time derivative. 

"(P ^ jj" After rescahng we may assume that (e*^)t>o is contractive. It is known that 
the invariance of y under {e^^)t>o is equivalent to invariance of y under R[X,A) for all 
A G IR large enough, see e.g. [26, Prop. 2.1]. On the other hand, [e'^*^)t£V is also a (unitary, 
hence contractive) Co-(semi)group that satisfies e'^^y cy ii and only if R{\ iA)y C y for 
A large enough, i. e., if and only \i iR {^ h , y G y for A large enough. In fact, the resolvent 
set of A contains an open sector of C which contains E, i.e., it contains the closed right 
half plane (with the possible exception of the origin). Then, for any Aq,/^ G S such that 

~ -^ol < [[^(-^o, ^)||~^ it is possible to develop the resolvent operator i?(/i. A) as a power 
series centered at Aq, i.e., 

OO 

i?(Ai,v4)-5](Ao-Ai)"i?(Ao,v4)"+i. 

Let now y be invariant under (e*^)t>o. Then y is invariant under R[\q,A) for some Ao, 
i. e., i?(Ao, A)y G y for all y & y. Since 3^ is a closed linear subspace, one obtains 

C30 

i?(Ai,v4)2/ = ^(Ao-Ai)"i?(Ao,v4)"+iyG3^ for all y G 3^ and |^-Ao| < \\R{\o,A)r\ 

and therefore y is invariant under i?(/i. A). This shows that the subset of the resolvent set 
for which i?(A, A)y C F is open. Moreover, it is relatively closed since y is closed. As a 
consequence, y is invariant under i?(/i. A) for all in the unbounded connected component 
of the spectrum containing Ao, and therefore also for all iA, A G IR large enough. By the 
representation of the semigroup in terms of the resolvent this shows that y is invariant under 
the unitary group (e'*'^)tgiR. 

"O lb])" This can be proved in the same spirit as the implication "(jb]) =4> Q". □ 
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Remark 5.4. Careful examination of the proof above shows that a(/, f) — a (e**^^/, e**''^/) 
and admissibihty of K are equivalent to 

(5.3) Vke*"^ = e*^>K for all t > 0, 

even if A is not self-adjoint. Definition [521 can thus be generalized by saying that Vk reflects 
a symmetry of a (possibly non- self- adjoint) parabolic network problem if ()5.3p holds. 

Thus, y is invariant under the action of (e**'*)t>o if and only if the associated orthogonal 
projection V is admissible and the orthogonality condition is satisfied by M, i. e., if and only 
if y is invariant under the action of (e*"^)t>o. In particular, for a star graph G Proposition 
14.71 (1) yields that there are nontrivial invariant subspaces of the above form if and only if 
there is a pair of edges with the same length. 



References 

[1] W. Arendt. Semigroups and evolution equations: Functional calculus, regularity and kernel estimates. 

In CM. Dafermos and E. Feireisl, editors. Handbook of Differential Equations: Evolutionary Equations 

- Vol. 1. North Holland, Amsterdam, 2004. 
[2] W. Arendt, C. J.K. Batty, M. Hieber, and F. Neubrander. Vector-valued Laplace Transforms and Cauchy 

Problems, volume 96 of Monographas in Mathematics. Birkhauser, Basel, 2001. 
[3] A. Arvanitaki. Effects evoked in an axon by the activity of a contiguous one. J. Neurophysiol, 5:89-108, 

1942. 

[4] S. Binczak, J.C. Eilbeck, and A.C. Scott. Ephaptic coupling of myelinated nerve fibers. Physica D, 
148:159-174, 2001. 

[5] H. Bokil, N. Laaris, K. Blinder, M. Ennis, and A. Keller. Ephaptic interactions in the mammalian 

olfactory system. J. Neurosci., 21:21:RC173, 1-5, 2001. 
[6] J. Bolte, S. Cardanobile, D. Mugnolo, and R. Nittka. Symmetries in quantum graphs. In W. Arendt 

and W. Schleich, editors. Evolution, Information, and Complexity. Wiley. 
[7] J. Boman and P. Kurasov. Symmetries of quantum graphs and the inverse scattering problems. Adv. 

Appl. Maths., 35:58-70, 2005. 
[8] V.I. Burenkov. Sobolev Spaces on Domains. Texte zur Mathematik 137. Teubner, Stuttgart, 1998. 
[9] S. Cardanobile and D. Mugnolo. Qualitative properties of coupled parabolic systems of evolution equa- 
tions. Ulmer Seminare iiber Differentialgleichungen und Funktionalanalysis, 11:85-103, 2006. 
[10] K.-J. Engel and R. Nagel. One-Parameter Semigroups for Linear Evolution Equations, volume 194 of 

Graduate texts in mathematics. Springer- Verlag, New York, 2000. 
[11] P. Exner and P. " Seba. Free quantum motion on a branching graph. Rep. Math. Phys., 28:7-26, 1989. 
[12] M. Kramar Fijavz, D. Mugnolo, and E. Sikolya. Variational and semigroup methods for waves and 

diffusion in networks. Appl. Math. Optim., 55:219-240, 2007. 
[13] B. Gutkin and U. Smilansky. Can one hear the shape of a graph? J. Phys. A, 34:6061-6068, 2001. 
[14] G. R. Holt. A critical reexamination of some assumptions and implications of cable theory in neurobi- 
ology. PhD thesis, California Institute of Technology, 1998. 
[15] G.R. Holt and C. Koch. Electrical interaction via the extracellular potential near cell bodies. J. Comp. 

Neurosci., 2:169-184, 1999. 
[16] B. Katz and O.H. Schmitt. Electric interaction between two adjacent nerve fibers. J. Physiol., 97:471— 
488, 1940. 

[17] V. Kostrykin and R. Schrader. Kirchhoff's rule for quantum wires. J. Phys. A, 32:595-630, 1999. 
[18] T. Kottos and U. Smilansky. Quantum chaos on graphs. Phys. Rev. Lett., 79:4794-4797, 1997. 
[19] P. Kuchment. Graph models of wave propagation in thin structures. Waves in Random Media, 12:1-24, 
2002. 

[20] G. Lumor. Connecting of local operators and evolution equations on networks. In F. Hirsch, editor. 
Potential Theory (Proc. Copenhagen 1979), pages 230—243, Berlin, 1980. Springer- Verlag. 

[21] G. Lumer. Espaces ramifies et diffusion sur le reseaux topologiques. CR. Acad. Sc. Paris, 291:627-630, 
1980. 

[22] G. Major, J. D. Evans, and J. J. Jack. Solutions for transients in arbitrarily branching cables: I. voltage 

recording with a somatic shunt. Biophys. J., 65:423-449, 1993. 
[23] D. Mugnolo. Gaussian estimates for a heat equation on a network. Math. Z., 2007. 

[24] D. Mugnolo. A variational approach to strongly damped wave equations. In W. Arendt et al., editor. 

Functional Analysis and Evolution Equations. Birkhauser, 2007. 
[25] S. Nicaise. Spectre des reseaux topologiques finis. Bull. Sci. Math., II. Ser., 111:401—413, 1987. 
[26] E. Ouhabaz. Analysis of Heat Equations on Domains, volume 30 of LMS Monograph Series. Princeton 

University Press, Princeton, 2004. 
[27] J.-P. Roth. Spectre du laplacien sur un graphe. C. R. Acad. Sci. Paris Ser. I Math., 296:793-795, 1984. 
[28] S. Severini and G. Tanner. Regular quantum graphs. J. Phys. A, 37:6675-6686, 2004. 



SYMMETRIES IN STRONGLY COUPLED NETWORK EQUATIONS 



25 



[29] J. von Below. A chajracteristic equation associated with an eigenvalue problem on C^-networks. Lin. 
Algebra Appl, 71:309-325, 1985. 

Institut fur Angewandte Analysis, Universitat Ulm, Helmholtzstrasse 18, D-89081 Ulm, Ger- 
many 

E-mail address: stefano.cardanobile9uni-ulm.de 
E-mail address: delio.mugnolo9uni-ulm.de 
E-mail address: robin.nittkaQuni-ulm.de 



